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A generalized inversion of block T + H matrix is obtained for the first time. 
In a particular case when T + H matrix is invertible, the method allows to 
obtain its inverse matrix without an additional condition of invertibility of the 
corresponding T — H matrix. 
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1. Introduction 

In many applications, e.g. digital signal processing, discrete inverse scatter¬ 
ing, linear prediction etc., Toeplitz-plus-Hankel (T + H) matrices need to be 
inverted. (For further applications see (Tj and references therein). 

Firstly the T + H matrix inversion problem has been solved in 0 where it 
was reduced to the inversion problem of the block Toeplitz matrix (the so-called 
mosaic matrix). The drawback of the method is that it does not work for any 
invertible T + H matrix since it requires also invertibility of the corresponding 
T — H matrix. This drawback appeared also in A, A- Later on the drawback 
was put out (see, e.g. E9. El), moreover, the inversion problem was solved for 
the block T + H matrix 0,0 

The generalized inversion of Toeplitz-plus-Hankel matrices is of great interest 
to, e.g., Pade-Chebyshev approximations. The generalized inversion for matrix 
A is meant to be the matrix such that AA A = A. Our goal is to obtain 
the generalized inversion of block T + H matrices with help of the well-known 
method of reducing the block-diagonal matrix formed from the T+H and T—H 
matrices to the mosaic matrix (as in the work 0). 

We will need the generalized inversion for the block Toeplitz matrix which 
has been already found in, e.g. 0. It is shown in the present paper that there 
is no need for T — H matrix to be inverted: if the T + H matrix is invertible 
than the obtained generalized inverse matrix proves to be its inverse matrix. 

The paper is organized as follows. At first the basic definitions and the main 
results of the paper 0 are given, then our main theorem is formulated and 
proved. This theorem is demonstrated with an example in the end of the paper. 

This work was supported by Russian Foundation for Basic Research (RFFI), 
grant N 04-04-96006. 

2. The basic definitions and notations 

We are going to find a generalized inversion for the block Toeplitz-plus- 
Hankel matrix 
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with CLj,bj G C pxq . 

We denote a" m (z) = a- m z~ m + ... + o 0 + ... + a n z n , &o +m (z) = b 0 + b\z + 
... + b n + m z n+m and introduce an auxilary matrix function 


A(z) 


« n 6J +m («- 1 ) z n ~ m a n i m {z~ 1 ) \ 
al m (z) z-™b% +m (z) J ■ 


Obviously, A(z) 


n 

£ A jZ j , with Aj G C 2px2q and 


j _ / b n —j o, n — rn —j 

Uy bj+m 


( 2 . 1 ) 


Thus, A(z) is the generating function for the sequence of matrices A - m ,..., Ao, 

• • -1 A n . 

Later on we will need a generalized inversion for the block Toeplitz matrix 


T a 
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A\ 


A. i 

Ao 


A—m \ 

A—m+l 


\ A n A n _i 


A 


n—m 




It has been already found in [5[. In order to use this result we should 
introduce the notations of the essential indices and polynomials of the sequence 
A — m ,..., Ao ,..., A n . 

We include the matrix Ta = To into the family of the block Toeplitz matrices 


T k = 


( A k 

Ak+i 


Ak -1 
Ak 


• ■ A— m ^ 

■ ■ A—rn+l 


—m < k < 


\ A n A n —i ... A n — m —k J 


The matrices Tk are of the same structure and it is reasonable that they 
should be examined together. 

We are interested in right kernels of Tk- For the sake of convenience let 
us pass from the spaces ker R Tfc to the isomorphic spaces TV” of generating 
polynomials. 

To do this we define the operator a a acting from the space of rational matrix 
functions R(z) = Y^=- n r i z ^i r o e C 2qxl to the space C 2pxl according to 


<t r {R{z)} = A -j r j- 

j=—n 


By A/)f, k = —m,..., n, we denote the space of vector polynomials R{z) = 
r j z ^ r 3 e C 2?xl , such that 

<j r { z~ l R(z )} =0, i = k, k + 1,..., n. 

K is evident to be isomorphic to ker R T)t. 

It is convenient to put NH m _ 1 = 0 and denote by 7Vjf +1 the 2 (n + m + 2 )q- 
dimensional space of all vector polynomials in z with formal degree n + m + 1. 


2 






Similarly, one may define spaces Aff which are isomorphic to ker L T)t. We 
denote ker^ A = {y \ yA = 0}. 

Let us put also a = dim Aff rn and to = dim A/j£. 

We will say that the sequence A- m ,...,A n is left (right) regular if a = 
0 (u = 0). Otherwise, the sequence is not regular and a(u>) is its left (right) 
defect. The sequence is called regular if a = u> = 0. It is evident that a < 
2 q, to < 2p for the nonzero sequence. 

We will denote d R = dim A/"*, A R = d® — d R _ x , k = —m, ...,n+ 1. 

As it is proved in 0, for any sequence A_ m ,..., A n the following inequalities 
hold: 

a = A R _ m < A R _ m+l < ... < A* < A* +1 = 2 (p + q)-u>. 

It means that there are 2 (p+q) — a—ui integers p a < p a +i < ... < p 2 (p+q)-a-ui 
satisfying equations 


A* 

= ... = Af 

—m 

Va+1 

A R 

= ... = A* 

1 

R 

A R 

^2(p + g)-a;+l 

~ • • • — ^n+l 


= a, 

= i, ( 2 . 2 ) 

= 2 (p + q)-u. 


If the ith row in (I is absent, we assume //.» = Pi+\. Let us put also = 
... = H a = -m - 1 if a ± 0 and y 2 ( P + q )-u 1 +i = • ■ • = T 2 ( P +q) if w ^ 0. 

Thus there is a set of 2 (p + q) integers, satisfying 12.2|) . for any sequence 
A - m ,..., A n . We will call these integers as indices of the sequence. 

Now we will define the right essential polynomials of the sequence. 

It follows from the definition of Af(f that Afjf and zAfff are the subspaces of 
A f£ +1 , k = — m — 1,..., n, moreover, A fjf H zAf R = Af R _ 1 . Then 


Afk+i — (Af R + zAf R ) © 7~t% +1 , 


where Tt R +1 is the complement of A/jf + zAfjf to the whole Af R +1 . 

Obviously, dim 7d” +1 = A£ +1 — A®. Hence dim 7d” +1 ^ 0 iff k = /Lt*. In this 
case dim7-f(( +1 is equal to the multiplicity ki of the index \Xi. 


Definition 2.1. If a ^ 0 then any column polynomials R\(z ),..., R a (z) form¬ 
ing the basis of Afff m will be called right essential polynomials of the sequence 
A_ m ,..., Aq, ...,A n . They correspond to the index pi = —m — 1 with the 
multiplicity a. 

Any vector polynomials Rj{z ),..., Rj+kj-i(z) forming the basis for Tt R j+1 
will be called right essential polynomials of the sequence A_ m ,..., Aq, ..., A n . 
They correspond to the index fij with the multiplicity kj, a + 1 < j < 2 (p+q)—to. 


Similarly, one may define the left essential polynomials. 

There are 2 (p + q) — to right and 2(p + q) — a left essential polynomials of 
the sequence A_ m ,... , A n . 

If a ^ 0 or to 0 then there is a lack of essential polynomials. But actually 
we can always complement the number of right (when p < q) or left (when p > q) 
essential polynomials to 2 (p + q). (The complement procedure was described in 

0 ). 

Henceforth, for definiteness sake, we will suppose that we have got the full 
set of 2 (p + q) right essential polinomials, i.e. either to = 0 or p < q. 
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The set of the left essential polynomials could always be recovered with the 
help of the so-called conformation procedure of the right and left essential poly¬ 
nomials. Let us describe how for the given set of the right essential polynomials 
Ri{z), ..., R2( p + q ){z), Rj(z) £ C 2 ?xl [z] one can construct the conforming left 
essential polynomials L\(z),..., L 2 ( p + q ){z), Lj(z) £ C lx 2 p [z]. 

We introduce the matrix of the right essential polynomials 


7Z(z) = (R 1 (z)...R 2(p+q) (z)) 


and find the matrix polynomial a_(z) from the next decomposition 
A(z)7Z(z) = a-(z)d(z) — z n+1 /3 + (z). 


Here d(z) = diag [z Ml ,..., z M2 < p +‘ >)], (3 + {z ) is the matrix polynomial in z, ck_(z) 
is the matrix polynomial in z _1 . Both of them have sizes 2p x 2{p + q). 

Let U-{z) be the matrix polynomial in z _1 such that: 


U-{z) 


K-(z) \ 

a-{z) ) ’ 


with Tl-(z) = z~ m ~ 1 TZ{z)d~ 1 (z). 

The matrix polynomial U-(z) is shown in jSj to be unimodular, i.e. its 
determinant is equal to a constant. 

We pick the 2 (p + q) x 2p block C(z) out UZ 1 (z): 


Uz\z) = ( * C{z) ). 


The matrix polynomial 


£(z) 


( Li(z) 

\ R 2 (p+q )(z) 


turns out to be the matrix of the conforming left essential polynomials. 

The case when a = 0 or p > q may be considered in a similar manner with 
help of the left essential polynomials. 

Now we may present the formula (5.13) from (H| for the generalized inverse 
of T a : 



( Ro ■ 

. 0 N 


( Co 

• \ 

T't _ 

1 A ~ 

V Rm ■ 

■ Rq ) 

" 

{ 0 . 

• Co ) 


Here IZj £ C 2gx2 ( p+g \ Cj £ C 2 ^ p+9 ^ x2p are the coefficients of the matrix 
polynomials 7?.(z), £(z), respectively, and Rj(z), Lj(z) are the conforming right 
and left essential polynomials of the sequence A_ m ,..., A 0 , ..., A n . The gener¬ 
alized inversion for matrix A is meant to be the matrix A' such that AZAA = A. 

The matrix n is constructed as follows. Let Ai,..., A r be the distinct es¬ 
sential indices of the sequence ... ,Ao,, A n and let Vi. ... , v r be their 
multiplicities (iq + ...-(- v r = 2 (p + q)). 
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Then 


n = 


( n 0 
n, 


n_i 

n 0 


n_„ \ 

n^n+l 


\ n m n?n—i • • • n m _ n J 

Here n fc = 0 for -n < k < to, k ^ -Ai,..., -A r , n_ Aj = \\sjSik\\ 2 ^=i\ 


1, i — v\ + • • • + Vj —i + 1,..., v\ + • • • + Vj, 
0, otherwise. 


The following partition of the right essential polynomials Rj (z) will be useful 
for the generalized inversion of the T + H matrix: 


Rj(z) = 


R){z) 

Rj(z) 


Here R l f 2 G C 9Xl [z]. In similar way we partition the left essential polynomials: 


L j( z ) = ( L )( z ) ) ’ 


with L)’ 2 G £ 1x p[z\. 

Then the matrix of these essential polynomials may be represented as: 

n z ) = ( ) . £(*) = ( £ 2 (z) ) , (2-4) 

with R}’ 2 (z) G C ?x2(p+9) , £}' 2 {z) G C 2 (p+9) x p. 

3. Generalized inversion 


In the section we will present our main result. 
Let us denote 



/ ni . 

.. 0 \ 

1 I 

■ 

■ ■ V-n \ 

II 

E? 

\ ni ■ 

... ^ 

II 

\ o • 

. . ) 


where Rj, (C 3 k ) are the coefficients of the polynomials R?(D). We also put 
Htz 2 = JT-R . 2 ) He i = T Cl J. 

Theorem 3.1. Generalized inverses oftheT+H andT—H matrices are found 
by the formulas: 


(T ±H)' = - 2 (T Ul ±H K2 ) H ( T C2 ±H Cl ). (3.1) 

If T ± H is invertible (one-sided invertible), then (T ± H)^ is its inverse 
(one-sided inverse) matrix. 
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Proof. 

Let us construct a generalized inversion to Ta = To according to formula 

(IP) . 

We are going to pass from block Toeplitz matrix Ta to the mosaic matrix 


bn 

• • • b n -\-m 

dn—m 

. . . CL n 

bn— 1 

• • • ^n+m—1 

Q-n—m—l 

• • • d n —1 

bo 

• • • bm 

d—m 

do 

a 0 

• • • Q—m 

bm 

bo 

ai 

• • • O'— ra+1 

bm -\-1 

... bi 


• • • d n —m 

bn-\-m 

• • • b n 


At first, according to the block structure of Ai EH we partition each 
block column Xj of the matrix Ta into two block columns Xj,Xj with sizes 
2 pin + 1) X q: 

X i = ( xj xj ) . 

Then permute new block columns in Ta and construct the matrix 




( x\ ••• xi 

xf ... 

xl ) = 


/ b n 
a 0 

bn-\-1 

a~ 1 

• • • b n -\-m 

• • • d—m 

dn—m 

bm 

dn— m+1 • • • 

bm — 1 • • • 

dn 

bo 

bn— 1 

ai 

bn 

a 0 

• • • b n -\-m—l 

• • • Oj — 771+1 

dn—m—1 

bm+1 

dn—m • • • 

bm 

dn— 1 

bi 



bo 

\ d n 

bi 

dn— 1 

• • • bm 

• • ■ O'n—m 

d—m 

bn-\-m 

d— m+1 • • • 
bn-\-m— 1 • • • 

do 

bn 


This matrix is evident to be obtained by multiplying Ta on a permutation 
matrix P 2 . Then we will do the analogous permutation with block rows in TAP 2 . 
As a result, we will get the matrix P 1 T 4 P 2 , where Pi is a permutation matrix. 
The matrix P 1 T 4 P 2 coincides with Ma ■ 


Ma = P\TaP2- 

Thus we have passed from the block Toeplitz matrix Ta to the mosaic matrix 

M a . 

Taking into account that for a permutation matrix P the equality P^ 1 = P* 
holds, we get the generalized inversion for Af 4 : 

M\ = P\T\Pt 

Let us specify the structure of factors in M\ = P|T^Pf. The operations 
which P 2 has done with the block columns of Ta, the matrix P.| now will carry 
out with the block rows of the matrix 

/ n 0 
\ n m 
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Thus 
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... 0 
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/ 

IZo 

' “I 
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nl 

... nl 




' : 


K 

... 0 


V 

7Zm ■ 

■ no y 








\ n 2 m 

... n'l 

J 


Tn, 

Tn 2 


where IZj’ are the coefficients of the matrix polynomials lZ 1,2 (z), presented in 

(IPl . 

( Cq ... C- r 


Similarly, we have 


/ n 


. p ‘ = 

\ 0 ... C 0 

rl r2 r2 

• A * / — n n 


= ( T Cl T C2 ) . 


Then 


V 0 ... Cl 0 ... Cl 

M a=(tI] ) n ( T ^ )■ 


Let us apply now the well-known method [2] of reducing the mosaic matrix 
Ma to the block-diagonal matrix formed from the Toeplitz-plus-Hankel and 
Toeplitz-minus-Hankel matrices: 


M a = 


1 ( J J 


2 V I -I 


T + H 0 
0 T-H 


We obtain that the matrix G = l 


I J 
-I J 




I J 
-7 J 

J J 
I -7 


l ( -rlYfri ) " (TcA + Tc, Tc,J-T i, ) 


is the generalized inversion for the matrix 
Let us present G in the block form 


T + H 0 
0 T-H 


G = 


G n G 12 
G 2 i G22 


with Gtj € C( m+1 )« x ( n+1 >. 
It is easy to get that 


Gn = - (T Wl + H U2 ) n ( T C2 +H Cl ), 


G 22 = 2 (Tn, -H n2 )H {T c , - H Cl ) 
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are generalized inverses to T + H (T — H). 

The theorem statement concerning the invertibility (one-sided invertibility) 
is evident. 

The theorem has been proved. 


Given T ± H matrices are block matrices with the sizes of their blocks p x q. 
The factors in the inverse formulas El have blocks with sizes q x 2(j> + q), 
2(p + q) x p. The compact form of the generalized inversion is in many respects 
because of such factors sizes. Sometimes it is convenient to have a formula for 
a generilized inversion where factors have blocks with sizes q x q,q x p,p x p. 
In order to obtain it we partition 1Z(z) = (i?i(z)... R 2 (j,+ q ) (z)) into blocks: 

-t-w , _ ( Till TZl2 R-13 ^14 \ 

{z) ~ V n 21 n 22 n 23 n 24 J • 


Here IZij have the sizes q x q for i. j = 1, 2, and q x p for i = 1, 2, j = 3,4. 
Then we will do the analogous partition with the matrix of the left essential 
polynomials: 


C{z) 


£11 

£12 ^ 

£21 

£22 

<£31 

£32 

£41 

£42 ) 


here Cij have the sizes q x p for i,j = 1, 2, and p x p for i = 3,4, j = 1, 2. 

Let us also partition D = diag[z Ml ... z^(p+<>>] = {d\ d 2 d 3 d ^), where d\- 2 
are diagonal matrices with the sizes q x q and ^ 3,4 are ones with the sizes p x p. 


We denote foi 

i,3 = 1, ■ 

..,4 





(K •• 

' 0 ^ 


( £0 ■ 

• £_„ \ 

T K« = 


•«->0 

... fV 

> T £ij = 

v 0 . 

.. £” / 


Then it is easy to see that 


(T±H)^= l - 


4 


j =1 


^2 H n 2j TrjHc n ± 


3 =1 


( 4 4 

X+ X lll <:>.,~i I c, 

3 =1 3 =1 

where we denote Tc,uJ = Hc,n and nj are the matrices constructed by dj with 
the same manner as n by d. 


4. An example 

Let us demonstrate the theorem with an example. We will find generalized 
inversions of the following T + H and T — H matrices: 


/ 

1 

0 

-1 

1 



/ 

1 

0 

-1 

1 

\ 


1 

1 

0 

-1 

± 


0 

-1 

1 

0 


1 

1 

1 

0 



-1 

1 

0 

0 


V 

-1 

1 

1 

1 

) 


V 

1 

0 

0 

1 

) 
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T±H = 












Our calculations show that the indices of the sequence A- 3 ,..., A 3 are equal 
to = —1, H 2 =0, H 3 = /X 4 = 1. The matrix IZ(z) of the right essential 

polynomials and the matrix C{z) of the conforming left essential polynomials 
are: 

n\z) 


k( z ) = ( n 2 {~z] ) ’ £( z ) = (£ 1 ( z ) £ 2 ( z ) ) 


where 


U 1 (z) = -(l 2 2 0 ) + ( —4 5 3 -l)z+ (2 1 0 0 ) z 2 - 

- ( 11 4 2 0 ) z 3 + ( 0 4 1 0 ) z 4 + ( 0 0 0 1 ) z 5 , 

K 2 {z) = ( 11 -4 —2 2) -( 2 11 5 -2 ) *+( 4 0 0 0 ) z 2 + 

+ ( 1 0 0 0 ) 2 3 + ( 0 0 1 0 )z 4 , 




( -4 \ 


( 4 \ 


f -1 \ 

1 


-12 

I 

0 

~-l 4- 

6 

180 


4 

1 

12 

Z 

-23 



^ 0 ) 


l 0 ) 


1 0 J 


z ~ 2 - 


( 21 \ 


( —18 \ 


( 16 \ 


-6 

z~ 3 4 - 

-84 

-" 4 4 - 

0 

z ~5 

35 

Z 1 

76 

^ 1 

0 


V 180 J 


l 0 j 


l 0 ) 



C 2 { Z ) = 


180 


( 25 \ 

30 

-25 

V 180 ) 


( 26 \ 

108 

-126 

V 0 ) 


z~ x + 


/-M 

30 
-43 

V 0 ) 


z~ 2 + 


/ -6 \ 


( 10 \ 


/ 4 \ 


-24 

4-3 

24 


0 

z~ 5 

32 

z — 

26 

z — 

0 

z 

V 0 j 


l 0 ) 


'v 0 J 



We construct the matrix 


n = 


( n 0 n_i n _ 2 n_3 ^ 

iii n 0 n_i n _ 2 

n 2 iii n 0 n_i 

V n 3 n 2 iii n 0 j 


where 


n n = 


( 0 

0 

0 

0 ^ 


f 0 

0 

0 

0 \ 

0 

1 

0 

0 


0 

0 

0 

0 

0 

0 

1 

0 

, n_i — 

0 

0 

0 

0 

V 0 

0 

0 

0 J 


^ 0 

0 

0 

1 ) 


n x = 


/ 1 0 0 0 \ 
0000 
0000 
V 0 0 0 0 / 

and the matrices II± 2 , n ±3 are zero matrices. 
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Then we construct the matrices , T-r 2 with help of the coefficients of 
1Z 1 , 1Z 2 and matrices , Tc 2 with help of the coefficients of C l , C 2 . 

We use the formulas ED for generalized inverses of the T ±H matrices and 
have 


(T + H ) f 


(T - H) ] = 




( 5 

10 

0 

° ^ 



1 


2 

-4 

12 

-4 



20 

-4 

8 

-4 

8 





l 1 

-2 

-4 

8 ) 




/ 

—113 

2 

50 

-32 

\ 



88 

8 

20 

52 



) 


-134 

-4 

80 

64 




V 

17 

-158 

10 

16 


) 


One may conclude that (T + H )' is the inverse matrix for the T + H matrix 
and (T — H)^ is the generalized inverse for T — H. 
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